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Abstract 

This research aimed to investigate the role of diverse transcription factors (TFs) and to delineate gene regulatory net- 
works directly in conifers at a relatively high-throughput level. The approach integrated sequence analyses, transcript 
profiling, and development of a conifer-specific activation assay. Transcript accumulation profiles of 102 TFs and 
potential target genes were clustered to identify groups of coordinately expressed genes. Several different patterns 
of transcript accumulation were observed by profiling in nine different organs and tissues: 27 genes were preferential 
to secondary xylem both in stems and roots, and other genes were preferential to phelloderm and periderm or were 
more ubiquitous. A robust system has been established as a screening approach to define which TFs have the ability 
to regulate a given promoter in planta. Trans-activation or repression effects were observed in 30% of TF-candidate 
gene promoter combinations. As a proof of concept, phylogenetic analysis and expression and trans-activation data 
were used to demonstrate that two spruce NAC-domain proteins most likely play key roles in secondary vascular 
growth as observed in other plant species. This study tested many TFs from diverse families in a conifer tree species, 
which broadens the knowledge of promoter-TF interactions in wood development and enables comparisons of gene 
regulatory networks found in angiosperms and gymnosperms. 

Key words: Conifer, expression pattern, Picea glauca, secondary cell wall, somatic embryogenesis, trans-activation assay, 
transcription factor, xylem. 



Introduction 

Plants produce wood through a process of secondary xylem 
development and growth with four major steps: initial cell divi- 
sion, cell expansion, cell-wall thickening, and programmed 
cell death (Plomion et ctl., 2001). Characteristic structures are 
thus formed that impart strength to support the plant body 
and the ability to transport water and nutrients efficiently 
from the roots to other organs (Miyashima et al, 2013). 
Secondary vascular growth represents a key adaptation and 



a major carbon sink in woody perennial plants; however, it is 
also observed on a smaller scale and studied in herbaceous 
plants. In Arabidopsis, it has been shown that the events of 
underlying primary and secondary vascular development are 
under tight regulatory control and involve transcription fac- 
tors (TFs) from the R2R3-MYB (MYB), basic helix-loop- 
helix, KNOTTED 1 -like homeobox (KNOX), homeodomain 
leucine zipper (HD-Zip), and NAM, ATAF1/2, and CUC2 
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(NAC)-domain families (Demura and Fukuda, 2007; Wang 
and Dixon, 2012). 

Recent studies have begun to unravel a transcriptional net- 
work regulating xylem development and growth. For example, 
the NAC-domain proteins are plant-specific TFs with several 
members playing a central regulatory role in xylogenesis, fibre 
development, and wood formation (Zhong et al, 2010a). 
Several NAC-domain genes are preferentially expressed in 
the cambial zone or xylem cells such as vascular elements and 
in interfascicular fibres (Demura and Fukuda, 2007; Zhong 
et al, 2010a; Wang et al, 2011; Wang and Dixon, 2012). 
In Arabidopsis, NAC-domain proteins (SND1, NSTs, and 
VNDs) orchestrate a transcriptional cascade involving MYB 
regulatory proteins, ultimately activating the transcription 
of enzymes involved in secondary cell wall (SCW) assembly. 
A similar cascade involving functional NAC-domain genes 
called WNDs (Wood-Associated NAC Domain TFs) and 
MYBs has been described in poplar and eucalyptus (Zhong 
and Ye, 2009; Zhong et al. , 201 1). 

The vascular tissues of angiosperms and gymnosperms 
have a number of anatomical and biochemical features in 
common and also many distinctive properties (Pallardy, 2008). 
These range from cell types and secondary cell macromol- 
ecules (both lignins and hemicelluloses) to intercellular con- 
nections and growth responses to stress (e.g. reaction wood). 
Conifers are a dominant group among the gymnosperms that 
may serve for comparative developmental studies in addition 
to being economically and ecologically important. Genomic 
initiatives in conifers have developed large gene sequence 
data sets, gene maps, and expression profiles (reviewed by 
Mackay et al, 2012) and more recently the complete genome 
sequences of two spruce species (Birol et al, 2013; Nystedt 
et al, 2013) and a pine (http://www.pinegenome.org/pineref- 
seq/). These information resources represent a foundation to 
accelerate the discovery of regulatory mechanisms of second- 
ary growth and wood formation. 

A class of transcriptional regulators that has been well stud- 
ied in conifer vascular development are the R2R3-MYBs. The 
loblolly pine {Finns taeda L.) genes PtMYBl and PtMYB4 
have been shown to be expressed in differentiating secondary 
xylem and to activate transcription through specific binding 
to AC motifs from promoters of most of the phenylpropa- 
noid pathway genes (Patzlaff et al, 2003a,b). When over- 
expressed in Arabidopsis and tobacco, these MYBs induced 
ectopic lignification (Patzlaff et al, 2003a; Newman et al, 
2004). Several other R2R3-MYB genes have been character- 
ized in white spruce (Picea glauca (Moench) Voss), loblolly 
pine (Bedon et al, 2007, 2010; Bomal etal, 2008) and in mar- 
itime pine (Pinus pinaster Ait.; Craven-Bartle et al, 2013). 
Sequence and phylogenetic analysis showed conservation in 
gene family structure between conifers and Arabidopsis. 

Research groups, including the current group, have over- 
expressed R2R3-MYBs and other TFs in transgenic conifer 
trees. Constitutive expression of PtMYBl or PtMYBS in 
spruce lines resulted in ectopic SCW deposition and increased 
lignin accumulation (Bomal et al, 2008). Craven-Bartle et al. 
(2013) independently showed that a Pinus pinaster homo- 
logue of MYB8 is a candidate regulator of phenylpropanoid 



metabolism and lignin synthesis genes. PtMYB8 is closely 
related to AtMYB61 and to AtMYB46, which are downstream 
of NAC-domain regulators (Zhong et al, 2010a), suggesting 
that PtMYB8 may be part of a transcriptional network con- 
trolling SCW deposition in conifers (Bomal et al, 2008). 

Nevertheless, the production of stably expressing conifer 
lines takes several months, and selection of lines with stable 
TF overexpression could be difficult to obtain due to phyto- 
toxic side effects or other pleiotropic effects. This group has 
developed a rapid transient transformation system for coni- 
fers to circumvent the problems related to the selection of 
specific appropriate transgenic lines (mostly due to integra- 
tion position effects) and nonspecific activation of genes or 
pathways through nontarget interactions. This article reports 
the development of a method for functional testing of TFs 
using a robust Agrobacterium-based transformation proto- 
col using embryonic spruce cells, which are easily cultured to 
facilitate the screening of a large number of genes directly in 
a conifer. This system enables the simultaneous introduction 
of TF-expression vectors and promoter-reporter gene con- 
structs into spruce and evaluation of resulting gene expression 
on a weekly basis. Over 600 functional assays were carried 
out with 59 full-length cDNAs encoding TFs cloned into 
plant expression vectors and putative promoter sequences 
for 12 candidate genes from Picea glauca. The portfolio of 
TFs and promoters was developed from the Picea glauca gene 
catalogue (Rigault et al, 201 1). The current work reports the 
identification of several TF-promoter combinations show- 
ing ?ra;M-activation using this experimental system. Sequence 
analyses, expression clustering profiling and clustering, and 
/raws-activation results were integrated to identify a putative 
spruce NAC TF proposed to control transcriptional activa- 
tion for genes involved in SCW formation in Picea glauca. 

Materials and methods 

Plant material 

For transcript profiling, 3-year-old white spruce seedlings obtained 
from open-pollinated seed lots were transferred to 2-1 pots and 
grown in a greenhouse under natural photoperiod from mid-April 
to the end of June. The plants were watered every 2 d on average and 
were fertilized weekly with 20 g T 1 N-P-K. Destructive sampling of 
48 plants involved four time points: 6 am and 3 pm on 29 June (day 
1) and 30 June (day 2), 2010. Nine vegetative tissue samples were 
collected from each tree: apex (AX, top 5 mm, main stem with nee- 
dles), secondary xylem (XS, main stem formed in the previous year), 
root xylem (XR, woody tissue of the largest roots), phelloderm (PS, 
bark and secondary phloem of the main stem formed in the pre- 
vious year), root periderm (PR, from the largest roots), elongating 
shoot (ES, annual growth without needles), root tips (RT), young 
needles (YNS, foliage on annual growth of the main stem; YNB, 
from branches). At each time point, 16 plants were sampled and the 
tissues from four trees were pooled to form a biological replicate 
(n = 4). Samples were quickly frozen in liquid nitrogen and stored at 
-80 °C until further use. 

RNA extraction, cDNA preparation, and quantitative PCR 
analysis 

Total RNA was isolated following the method of Chang et al. (1993) 
with modifications described previously (Pavy et al, 2008). RNA 
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concentration and integrity were determined using a NanoDrop 
1000 (Thermo Scientific, Wilmington, DE, USA) and an Agilent 
2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA), 
respectively. Complementary DNA (cDNA) was synthesized 
with the Superscript First-Strand Synthesis System for RT-PCR 
(Invitrogen, Carlsbad, CA, USA), following the manufacturer's 
instructions with minor modifications. Briefly, 1 ug total RNA was 
reverse transcribed using an anchored oligo d(T), and a GFP spike- 
in was added as an internal control of the reverse transcription. The 
cDNAs were diluted 1 :4 in RNase-free nanopure water. Quantitative 
PCR was essentially as per Boyle et al. (2009) with some modifica- 
tions. Briefly, PCR mixtures contained either QuantiFast SYBR 
Green PCR kit (Qiagen, Germantown, MD, USA) or LightCycler 
480 SYBR Green I Master (Roche, Basel, Switzerland) and were 
composed of 1 x master mix, 300 nM of each gene-specific primer, 
and 5 ul diluted cDNA in a final volume of 15 ul. Gene-specific 
primers (Supplementary Table SI available at JXB online) were 
designed with Primer3Plus software (http://www.bioinformatics.nl/ 
cgi-bin/primer3plus/primer3plus.cgi), verified for self-complemen- 
tarity with the Oligonucleotide Properties Calculator (http://www. 
basic.northwestern.edu/biotools/oligocalc.html), and for specific- 
ity against the Picea glauca gene catalogue (Rigault et al, 2011). 
Reaction specificity was verified by the presence of a single ampli- 
fication product based on melting curve analysis. The number of 
molecules (transcripts) were calculated for each sample using linear 
regression of efficiency (Rutledge and Stewart, 2008) adapted for 
Excel (Boyle et al, 2009), normalized using a ratio based on the 
geometric mean of six Picea glauca reference genes and the GFP 
spike-in (HM1 5 1400.1) and transformed to a log 2 scale for statistical 
analysis. Melting curve analysis was performed to ensure the ampli- 
fication of a single product. Samples whose amplification of primer- 
dimer or nonspecific amplification was equivalent to 50% or more of 
the amplification of the gene were removed from the analysis. 

The reference genes were: elongation factor la (EFl-a, BT102965), 
cell division cycle 2 (CDC2, BT106071), ribosomal protein L3A 
(BT1 15036), eukaryotic initiation factor 4E (BT1 12014), ubiquitin- 
conjugating enzyme (BT 109864), and core histone H2AIH2BIH3I 
H4 (BT1 16867) (Supplementary Table S2). The stability and the 
variation of the reference genes were analysed with the geNorm 3.5 
algorithm (Vandesompele et al, 2002) before using them for nor- 
malization (Supplementary Table S7). 



Analysis of expression data 

Normalized quantitative PCR data of the 102 different genes were 
obtained for day 1 and were analysed to compare the morning and 
afternoon samples (for each gene and tissue separately) by using 
Student's t-test, with Bonferroni's step-down correction for multi- 
ple testing (Holm, 1979). Further analyses only included genes with 
no significant difference between morning and afternoon samples 
(adjusted P-value < 0.05), which were combined (within a tissue) for 
statistical tests. The effect of the tissue type was tested by using an 
individual one-way analysis of variance for each gene; P-values were 
adjusted with Bonferroni's correction. Duncan's multiple range test 
was used to determine preferential expression among the nine tissues. 
The analyses were conducted in the R statistical environment release 
2.12.1. In addition, 15 genes were tested for consistency between 
day 1 and day 2. Strong and statistically significant Pearson's coef- 
ficients (mean 0.86) were obtained between days for the 13 genes that 
varied between tissues, and no interactions were observed between 
day and tissue effects, even though the day effect was significant for 
some genes (Supplementary Table S3). In contrast, the two genes 
that did not vary between tissues were found to vary significantly 
between days and gave low (nonsignificant) Pearson's coefficients. 
This work analysed only samples from day 1 given that tissue dif- 
ferential expression was repeatable over this time period. 

The genes with significant tissue effect were clustered by using 
complementary methods in the MultiExperiment Viewer soft- 
ware release 4.2 (http://www.tm4.org; Saeed et al, 2003). Prior to 



clustering, the scale of the data was transformed from numbers of 
molecules to relative expression values: for each gene, the tissue with 
the most RNA molecules was set to 100, the threshold of detection 
was set to 0, and a relative value (between 0 and 100) was calcu- 
lated for other tissues based on their ratio to the highest sample. 
Hierarchical clustering of genes and tissues was used to represent 
tissue-preferential transcript accumulation. The self-organizing tree 
algorithm (SOTA) was used to cluster genes and define expression 
profiles according to the following parameters: maximum cell diver- 
sity 0.01, Euclidean distance, and a maximum of 12 cycles. Each 
gene was statistically tested for fit against the 12 SOTA vectors 
by using Pavlidis template matching (PTM) analysis with default 
parameters and a P-value of 0.05. 

Isolation of 5' -genomic sequences and construction of vectors 

The 5'-upstream sequence of each candidate target gene's cod- 
ing region was identified using the Universal GenomeWalker kit 
(Clontech, Mountain View, CA, USA) as previously described 
(Germain et al, 2012). Amplified fragments were cloned using the 
TA cloning kit (Invitrogen) and electroporated into Escherichia 
coli XL 1 -Blue. Promoters were subsequently amplified from Picea 
glauca genomic DNA using gene-specific primer pairs designed from 
the fragment obtained by genome walking (Supplementary Table 
S2) and then cloned and sequenced. The putative promoter regions 
and 5'-UTR of each candidate gene were submitted to GenBank 
(Table 2). The digested fragments (Xbal-BamHT) were inserted into a 
modified pMJM vector (Levee et al. , 2009) to create promoter::GUS 
fusions as previously described in Bedon et al, (2009). These 
fusion fragments were inserted into a modified pCAMBIA2300 at 
the &<?8387I site where the hygromycin resistance gene is removed 
and replaced by the silencing inhibitor pl9 protein-encoding gene 
(Voinnet et al, 2003) to improve transgene expression efficiency. 
TF-expression vectors were obtained by PCR amplification using 
cDNA full-length coding sequences. PCR products were transferred 
to a pCAMBIA2300 expression vector using the Gateway system 
(Invitrogen), where they were flanked by the ubiquitin promoter and 
the 35S terminator. The reporter and expression vectors were trans- 
ferred to Agrobacterium tumefaciens AGL1. 

Transient trans-activation transformation and quantification 
analysis 

A rraw^-activation assays system based on Agrobacterium cotrans- 
formation using Picea glauca somatic embryogenic cells was devel- 
oped using an approach similar to a previously described one for 
Arabidopsis (Berger et al , 2007). This transient transformation assay 
is derived from the initial steps of the stable transformation proce- 
dure described for spruce (Klimaszewska et al. , 2004) with conditions 
optimized to maximize and maintain the level of transient expression 
observed up to 6 d. Briefly, a liquid culture of Picea glauca embryo- 
genic cells (line PG653) was maintained and subcultured weekly. 
Following an overnight incubation, cells from the two Agrobacterium 
cultures (promoter-reporter with TF or empty expression vectors) 
were harvested by centrifugation, resuspended in embryogenic liquid 
culture media, and each added to 10ml of 1-d-old suspension culture 
to a final OD of 0.3. The cocultures were then placed in a shaker 
for lh at 21 °C after which the cells were evenly vacuum filtered 
onto four separate 4.3-cm filter papers. Each filter was then placed 
in a separate Petri dish on two 18-cm filter papers wetted with 7 ml 
liquid media supplemented with 100 uM acetosyringone. The Petri 
dishes were then sealed and kept in the dark at 21 °C After 6 d of 
coculture, each filter paper was vacuum pulsed to remove excess liq- 
uid and a histochemical GUS assay was performed (Klimaszewska 
et al. , 2004). The level of /raws-activation on each filter was visually 
evaluated compared to the control without TF. For each of the TF 
promoter combinations, the observed reporter gene expression was 
rated based on the number of filters showing increased GUS blue 
staining. Thms-activations were rated as positive when visually at 
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least three of the four filters showed clearly increased expression of 
GUS by histological staining (relative to the control), or likely posi- 
tive when only two filters showed clearly increased GUS expression. 
In the absence of increased GUS expression on any of the four filters, 
or when only one filter was positive, the test was evaluated as show- 
ing no rra/K-activation. For those promoters alone, which showed a 
sufficiently high level of background GUS expression, promoter-TF 
interactions that caused a reduction in observed promoter activity 
were evaluated as causing a downregulation. 

Sequence comparison and phylogenetic analysis 

Full-length nucleotide sequences for 27 NAC genes were obtained 
from this group's Picea glauca gene catalogue (Rigault et al., 201 1). 
Class IIB NAC-domain subfamily was defined using Arabidopsis 
and poplar members (Ohtani et al, 201 1). The NAM/NAC domain 
from deduced amino acid sequences was used for multiple alignment 
using ClustalW with MEGA 4.0.2 software (Tamura et al, 2007), 
as previously described (Bedon et al, 2010). Two neighbour-joining 
trees were constructed using MEGA 4.0.2 software using the fol- 
lowing parameters: Jones-Taylor-Thornton substitution model, 
gamma parameter of 1 .0, and pairwise gap deletion. Bootstrap con- 
sensus trees were deduced from 1000 replicates. 



Results 

Transcript profiling and expression cluster analysis of a large 
number of TFs and their potential target genes allowed this 
work to identify groups of coordinately expressed genes with 
respect to wood formation. A portfolio of gene constructs of 
full-length coding sequences for several classes of TFs was 
developed as well as the promoter sequences of selected genes 
involved in different physiological processes related to wood 
formation. TF/promoter interactions were identified using 
a transient transformation system developed in an embryo- 
genic white spruce cell culture. 

Clustering of transcript profiling data 

High-quality transcript accumulation data were obtained for a 
set of 102 TFs and candidate genes in nine different vegetative 
tissues or organs from Picea glauca (Supplementary Table S6). 
A large majority of the genes tested (82/102, 79%) showed sta- 
tistically significant variations between the nine different sam- 
ple types (ANOVA P<0.05 with correction for multiple tests). 
A much smaller proportion (19%) varied significantly between 
the morning and afternoon in one or more sample types 
(Student's t-test, P<0.05 with correction). There was little over- 
lap between the tissue preferential genes and the genes that var- 
ied between the morning and afternoon samples. Hierarchical 
clustering was applied to the 82 genes that varied between tissue 
and organ samples to visualize putative groups of coexpressed 
genes. The genes were clustered into two major expression 
groups (Fig. 1, left panel): one cluster contained genes with a 
clear preferential accumulation either in the secondary xylem 
of the stems and roots or in the phelloderm and periderm 
(Fig. 1, left panel, bottom). A second and larger cluster of 
genes included a variety of profiles, with moderately preferen- 
tial expression in various tissue types (Fig. 1, left panel, top). 

The expression data were analysed with SOTA followed by 
PTM to define statistically robust coexpression groups for the 



82 tissue preferential genes. The SOTA analysis grouped the 
genes into five major expression clusters comprising between 
six and 25 members (clusters 1, 2, 3, 7, and 9; Supplementary 
Fig. SI), and seven minor clusters containing between one and 
three members. The PTM analysis showed similarity among 
the SOTA clusters: 69 genes fitted at least one and up to five 
profiles with statistical confidence (P<0.05). Most notably, 27 
genes matched three or four of the clusters 9, 10, 11, and 12 
(although some also matched up to two other clusters). The 
TFs in this set of transcripts included eight of the R2R3-MYBs 
genes (PgMYBl, PgMYB4, PgMYBS, PgMYBl, PgMYB3 
PgMYB16, PgMYB19, and PgMYB27), two HD-ZIP III 
(PgHB7 and PgHBS), and three NAC (PgNAC-3, PgNAC-7, 
highly preferential to xylem, and PgNAC-8) that were pref- 
erential to secondary xylem tissues from stems and roots. 
They were coexpressed with candidate genes encoding lignin 
biosynthesis enzymes (Pg4CL, PgCAD, PgCCR-1, PgHCT, 
PgOMT-1) and cell-wall proteins (PgKORl, PgAGP-1). 

Development of a portfolio of conifer TFs and specific 
candidate gene promoters 

Based on predicted gene functions and expression data clus- 
tering, a portfolio of spruce TF-expression vectors and of 
promoter::GUS constructs for candidate genes was devel- 
oped for screening in functional assays. This work focussed 
on 13 TF families (Table 1) from the collection of full- 
length Picea glauca cDNAs in which 927 unique transcript 
sequences have been assigned to 34 plant TF families (Rigault 
et al. , 201 1). The complete list and annotations are presented 
in Supplementary Tables S4 and S5. Members of the selected 
families have been linked to development, including vascular 
growth and differentiation such as NAC, MADS, AUXMAA, 
to stress responses such as AP2, WRKY, or to both of these 
biological processes, such as b-Zip and MYBs. Transcript 
profiling and activity screening experiments covered an over- 
lapping set of 59 TF cDNA clones, each containing a com- 
plete coding sequence as determined by large-scale cDNA 
sequencing (Rigault et al. 2011) and targeted TF analyses in 
conifers (Bedon et al, 2007, 2010; Bomal et al, 2008). 

The choice of candidate genes for promoter isolation 
was based on functional data from previous experiments, 
expression data, and scientific literature indicating they may 
be downstream targets for the TFs under study. Overall, 
this study selected genes that have been linked to vascular 
growth and differentiation in trees, including those encoding 
for lignin biosynthesis enzymes (Pg4CL, PgCAD), proteins 
involved in cell-wall synthesis and remodelling (PgCesA-3, 
PgXTH8-l, PgTUA-1, alpha-tubulin-1), and transcriptional 
regulators (PgHB4, PgLIM-1, PgMYBl, PgMYB8; Table 2) 
based on several previous reports (Plomion et al, 2001; 
Peter and Neale, 2004; Bomal et al, 2008; Paiva et al, 2008; 
Bedon et al, 2010). The list also included wound-inducible 
and defence-response genes (PgCAD, PgDHS2; Bedon et al, 
2009). The putative promoters of these 12 genes were identi- 
fied as the region immediately upstream (5') of their coding 
sequences (Table 2). The lengths of the isolated promoters 
were from 931 to 3620 bp. 
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Fig. 1. Hierarchical clustering showing tissue-preferential transcript accumulation. The clustering was based on relative transcript accumulation values 
(0-100) determined by quantitative PCR with gene-specific primers and based on a mean number of molecules normalized for each tissue. Genes that 
varied between sampling time points (morning and afternoon) are not shown (grey boxes). XS, secondary xylem from the main stem; XR, secondary 
xylem from largest roots; ES, elongating shoot (without needles); PS phelloderm (bark and secondary phloem) of the main stem; PR, periderm from 
largest roots; AX, shoot apex; RT, root tips; YNS, young needles from the main stem; YNB, young needles from branches. 
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Table 1 . Selected TF families used for trans-activation assays 

AP2: APETALA2 domain; Aux/IAA: auxin/indole-3-acetic acid proteins; 
bZIP: basic region/leucine-zipper motif; GRAS: GIBBERELLIN ACID 
INSENSITIVE (GAI), REPRESSOR of GA1 (RGA), and SCARECROW 
(SCR); Homeobox: homeodomain; KNOX2: KNOX2 domain proteins; 
LIM: LIM domain; MADS: MCM1 , Agamous, Deficiens, Srf; MYB: 
Myb-like DNA-binding domain; NAC: NAM, ATAF1/2, and CUC2 
domain; WRKY: WRKY DNA-binding domain. 



Family cDNA Full-length Tested 

identified coding sequence interactions 



AP2 


63 


21 


93 


AUX/IAA 


26 


7 


24 


B-box zinc finger 


6 


6 


8 


bZIP 


31 


11 


64 


GRAS 


17 


4 


20 


Homeobox 


54 


9 


8 


KNOX2 


7 


3 


24 


LIM 


5 


4 


32 


MADS 


45 


18 


70 


MYB* 


122 


42 


201 


NAC domain 


36 


19 


84 


NASC 


7 


6 


44 


WRKY 


25 


12 


36 


Total 


444 


162 


708 



Table 2. Candidate genes used for promoter isolation 



Functional screening for TF-promoter trans-activation 

An Agrobacterium transient transformation method was 
developed for spruce and was used to test 664 combinations 
comprising a TF and a candidate target gene promoter. As 
illustrated in Fig. 2, the method used Picea glauca embryogenic 
cell cultures that were cotransformed with Agrobacterium: one 
culture containing a promoter: :GUS construct and the other 
culture either with a full-length coding sequence of a selected 
TF (to test for interaction) or an empty vector (control). 

Fig. 3 presents an assessment of GUS staining following 
each of the 664 cotransformation experiments, with each 
interaction tested in quadruplicate. Positive combinations 
are those in which the addition of the TF modified the level 
of GUS staining with a given promoter:: GUS construct. In 
total, this study observed 192 clear positive combinations, 
representing 29% of all the combinations tested. Promoter- 
TF combinations with positive GUS staining in two of the 
replicates suggested weak or conditional interactions and 
represented 13% of the tested interactions. No rra«^-activa- 
tion was observed when promoters were cotransformed with 
empty vectors (no TF), except for PgDHS2 and PgGASA5-l 
promoters, which showed significant and consistent back- 
ground GUS expression. A clear repressive effect (1% of 
interactions) by a few TFs was observed with the PgDHS2 
and PgGASA5-l promoters (Fig. 3). 



Size includes the promoter region and the 5'-UTR. 



Candidate 
gene 



Description 



Function 



Size (n) 



GenBank 
accession no. 



References 



PgDHS2 



Pg4CL 
PgCAD 

PgCesA-3 
PgHB4 

PgLIM-1 



PgMYBI 

PgMYB8 

PgGASA5-1 

PgSABATH2 

PgTUA-1 
PgXTH8-1 



3- Deoxy-7- 
phosphoheptulonate 
synthase 

4- Coumarate CoA ligase 
Cinnamyl alcohol 
dehydrogenase 
Cellulose synthase A 
Class III homeodomain 
leucine zipper (HD-Zip III 
LIM (UN11/ISL-1/MEC-3) 
TF family 

MYB domain TF 

MYB domain TF 

Gibberellic acid-stimulated 

Arabidopsis 

SABATH family methyl 

transferase 

a-Tubulin 

Xyloglucan 

endotransglycosylase/ 
hydrolase 



Enzyme of the shikimate pathway; provides 1 353 

precursors for monolignol and flavonoid metabolisms 



Lignin synthesis 
Lignin synthesis 

Cellulose synthesis in secondary cell walls 
Transcriptional regulation; involvement in primary and 
secondary vascular tissue pattern formation 
Transcriptional regulation of genes related to 
secondary cell walls 

Transcriptional regulation of genes related to 
secondary cell walls 

Transcriptional regulation of genes related to 
secondary cell walls 

Gibberellin-responsive protein; putative role in cell 
expansion in wood-forming tissue of Pinus pinaster 
Putative role in phytohormone methylation 

Cytoskeleton and cell-wall organization 
Modification of the xyloglucan-cellulose framework 
of plant cell walls; control of cell-wall expansion and 
strength 



1885 
1519 



2085 
2407 



1780 



995 
2179 



944 
931 



JN828804 



JN828803 
FJ428229 

KF824520 
KF834195 

KF834197 



3620 KF834198 
2332 a KF834199 



KF8341 96 
KF834200 



KF834201 
KF834202 



Bomal ef al. (201 4) 



Bomal ef al. (201 4) 
Bedon ef al. (2009) 

Taylor (2008)" 
Coteef al. (2010) 

Kawaoka and Ebinuma 
(2001)"; Demuraand 
Fukuda (2007) 6 
Bedon ef al. (2007); Bomal 
ef al. (2008) 

Bedon ef al. (2007); Bomal 
ef al. (2008) 
Paiva ef al. (2008)" 

Zhao ef al. (2009) 

Lloyd and Chan (2008)" 
Mellerowicz ef al. (2008)" 



includes one intron of 278 bp. 

b General references about putative gene function. 



(A) 
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Reporters 





1 


1 


p 


p19 


35S pro 




Candidate gene promoter 


GUS reporter 


NOSter 



Effectors 



Control (empty vector) 



Ubi pro 



(B) 



TF expression vector 



Ubi pro H Transcription factor 





Fig. 2. Overview of the components of the fra/is-activation method. (A) Schematic diagrams of the effector and reporter plasmids. The reporter 
constructs consisted of a GUS reporter gene driven Py the tested candidate gene promoter; p19 driven Py the CaMV 35S promoter was inserted into the 
same plasmid. The effector constructs consisted of the Gateway-inserted complete coding sequence of each TF driven by the maize ubiquitin promoter. 
The control (empty vector) consisted of the expression vector without TF insertion. (B) Representative positive frans-activation. Histochemical GUS assay 
was performed after 6 d of emPryogenic cells/agroPacteria coculture. Left: control assay of white spruce cells transformed with the PgMYB8 promoter- 
reporter construct and empty vector effector construct; right: white spruce cells transformed with the PgMYB8 promoter-reporter construct and the 
PgNAM- 7 effector construct. 



Gene expression levels of the TFs in spruce embryo- 
genie cells were obtained from the PiceaGeneExpress data- 
base (Raherison et al., 2012) and indicated that most of the 
TFs were among the 1 5% genes not detected by microarray 
analysis and only seven genes of the 33 tested were shown 
to be strongly expressed (Fig. 3, expression class >7). These 
observations indicated that the potential contribution of 
endogenous TF expression to fraws-activation was small 
overall considering that recombinant TFs are expected to be 
expressed at much higher levels than endogenous factors in 
transiently transformed cells. Overall, the data strongly sug- 
gest that the method enabled the identification of TFs that 
acted on the candidate gene promoters. 

Identification in conifers of a putative transcriptional 
regulation network similar to the SND1/VND cascade 

The development of a rapid method for screening TF-pro- 
moter interactions in a gymnosperm allowed us to survey the 
potential role of several members of the NAC-domain TF 
family in vascular growth in conifers. Functional studies of 
Arabidopsis NAC-domain proteins such as NST members 
(NAC SCW-thickening promoting factors) and SND mem- 
bers (SCW-associated NAC-Domain) have indicated that 



they are key transcriptional switches governing SCW biosyn- 
thesis and formation of vascular vessels, respectively (Demura 
and Fukuda, 2007; Yamaguchi et al. , 201 1 ; Wang and Dixon, 
2012). The current study investigated eight genes from the 
NAC-domain family to assess whether specific members 
with a predominant role in SCW formation may be found in 
gymnosperms. PgNAC-4 and PgNAC-7 gave positive inter- 
actions with the promoters of genes encoding enzymes of 
biosynthetic pathways related to SCW deposition (PgCAD, 
Pg4CL, PgDHS2, and PgXTH8-l) and genes encoding TFs 
expressed in wood-forming tissue (PgLIM-1 and PgHB4 
for PgNAC-4; PgMYBl and PgMYB8 for PgNAC-7). This 
observation is significant because transgenic Picea glauca 
plants overexpressing their pine orthologues, PtMYBl and 
PtMYB8, had enhanced lignification and developed ectopic 
SCWs, to a greater extent in PtMYB8 lines (Bomal et al, 
2008). The pine genes PtMYBl and PtMYB8 were also tested 
in this transient system and were able to positively regulate 
promoters from genes directly related to SCW synthesis such 
as PgCesA-3, Pg4CL, and PgXTH8-l (Fig. 3). PtMYBS was 
able to positively regulate the PgMYBl promoter, suggest- 
ing that PgMYBl may be downstream of PgMYB8. These 
results suggest that PgNAC-7 may be a TF that acts as an 
upstream regulator targeting PgMYBl and PgMYB8, which 
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^(Positive trans-activation 
| Ho trans-activation 
| Not tested 



Likely positive trans-activation 
[ Down-regulation 



ND, not detectable 
NT, not tested 



Fig. 3. Transient frans-activation assays leading to the identification of specific interactions. Each box represents the observation of the histochemical 
test compared with the control without TFs. The PgDHS2 and PgGASA5-1 promoters alone showed faint GUS staining and loss of expression (relative 
to the control) was observed with specific TFs (yellow). Relative gene expression levels for each TF based on PiceaGenExpress are indicated on the right. 
This database was developed from transcript profiles obtained for eight different tissue types coming from five independent experiments (Raherison et al. , 
201 2). Briefly, the genes were ranked based on their signal intensities within a tissue type and equally divided into 1 0 separate classes according to their 
signal intensity: class 1 : 1 0% lowest signal intensities; class 1 0: 1 0% highest signal intensities (see Raherison et al., 201 2 for details). None of the genes 
tested were ranked in class 10. 
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in turn positively regulate genes encoding enzymes of SCW 
formation. 

Identification of two conifer NAC-domain TFs similar to 
their orthologous proteins from angiosperms 

Previous studies reported that NAC-domain TFs involved in 
vascular development form the class IIB clade that includes 
NST1/2/3, VND1-7, SMB, and BRN1/2 (Kubo et al, 2005; 
Mitsuda et al. , 2005 ; Zhong et al. , 2006; Bennett et al. , 201 0). 
Here, a multiple alignment of NAC-domain sequences 



(N-terminus DNA-binding domain) and construction of 
a phylogenetic tree showed that PgNAC-4 and PgNAC-7 
are the only Picea glauca NAC-domain TFs (within the 
31 identified) that belong to the class IIB clade (Fig. 4A). 
PgNAC-4 was part of the BRN/SMB subgroup whose pro- 
teins were shown to play a role in root cap differentiation 
(Bennett et al, 2010). In contrast, PgNAC-7 showed higher 
homology with the VND family, particularly with VND4 
and VND5 (Fig. 4B). Members of the IIB NAC-domain 
family contain two characteristic motifs within their tran- 
scriptional activation domain located at the C terminus (i.e. 



(A) 



(B) 
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Fig. 4. Identification of conifer NAC-domain proteins similar to the vascular proteins from angiosperms. (A) Unrooted neighbour-joining phylogenetic 
tree of NAC-domain proteins identified in Picea glauca. Arabidopsis proteins were chosen as landmarks representing the vascular NAC-domain IIB 
subfamily (grey box). All Picea glauca NAC-domain proteins are indicated by open diamonds and those identified in the IIB subfamily are indicated by 
filled diamonds. (B) Unrooted neighbour-joining phylogenetic tree of NAC-domain proteins from Arabidopsis, poplar, and Picea glauca identified in the 
IIB subfamily. Picea glauca proteins are indicated by filled diamonds. (C) Multiple sequence alignment of the C-terminus domains between PgNAM-7 
and members of the VND subfamily; LP and WQ motifs are indicated. (D) Multiple sequence alignment of the C-terminus domains between PgNAM-4 
and members of the SMB/BRN subfamily; LP- and WQ-like motifs are indicated. At: A. thaliana; Pg: Picea glauca; Pt and Ptr: Populus trichocarpa 
sequences, published by Ohtani et al. (201 1 ) and Zhong et al. (201 0b), respectively. 
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the LP and WQ boxes necessary for their activity; Ko et al, 
2007). The protein sequence of the PgNAC-7 C termini 
contained canonical LP and WQ motifs (Fig. 4C) identical 
to those in AtVND4 and PtVNS03/PtrWND4A. PgNAC-4 
contained sequences related to LP and WQ motifs, mostly 
represented by the canonical proline and tryptophan as 
in the Arabidopsis genes AtSMB, AtBRNl, and AtBRN2 
(Fig. 4D). 

Expression profiling showed 8-120-fold higher transcript 
accumulation for PgNAC-7 in secondary xylem of both 
roots and stems compared to the other vegetative tissues and 
organs tested (Fig. 5). This finding is consistent with obser- 
vations reported for the VND4/5/6 subfamily in Arabidopsis 
(Kubo et al, 2005). Similarly, PgNAC-4 transcript levels 
were relatively uniform in most tissues but were 8-10-fold 
higher in the root tip (Fig. 5), also consistent with AtBRNl 
and AtBRN2 expression patterns in Arabidopsis (Bennett 
et al, 2010). Analysis of transcript accumulation for spruce 
TFs and candidate genes in the present study indicated that 
PgNAC-7 shared tissue profiles also found for genes directly 
involved in secondary growth. 



n 14 




Fig. 5. Transcript profiles of PgNAC-4 (top) and PgNAC7 (bottom) in 
white spruce. Transcript levels were determined by quantitative PCR with 
gene-specific primers. Data are mean standard deviation from three or 
four biological replicates, each comprised by pooling tissues from four 
trees, and are presented as log 2 number of RNA transcript molecules 
per ng of total RNA; differences of 2, 3, 4, or 5, represent 4-, 8-, 1 6-, or 
32-fold differences, respectively. ES, elongating shoot (without needles); 
PS, phelloderm (bark and secondary phloem) of the main stem; PR, 
periderm from largest roots; XS, secondary xylem from the main stem; XR, 
secondary xylem from largest roots; AX, shoot apex; RT, root tips; YNS, 
young needles from the main stem; YNB, young needles from branches. 



Discussion 

Coordinated expression of TFs and candidate genes in 
specific trans-activation patterns 

Despite the rapid development of genomic tools for coni- 
fers such as genetic maps, expressed sequence tag databases, 
and DNA microarray systems, functional characterization 
of candidate genes is challenging. Screening for mutants is 
an extremely tedious task and the long breeding time makes 
forward genetic approaches impossible. Populus has been a 
model tree species for studies of gene function because of ease 
of in vitro cultivation and genetic transformation, including 
activation tagging (Harrison et al. , 2007). For conifers such as 
pines and spruces, efficient genetic transformation procedures 
have been developed (Tang et al, 2001; Klimaszewska et al, 
2004). Nevertheless, many conifer genes have been studied in 
model plants such as Arabidopsis and tobacco because of the 
long regeneration times and the technical complexity of pro- 
ducing transgenic conifers (Patzlaff et al, 2003b; Newman 
et al, 2004). 

Only a limited number of genes have been functionally 
analysed by using a homologous expression system in conif- 
erous species. Some involve misexpression analyses of a gene 
encoding an enzyme to reveal its potential function in a meta- 
bolic pathway such as lignin (Wadenback et al, 2008; Wagner 
et al, 2009) and terpene biosynthesis (Hamberger et al, 
2011). Molecular mechanisms underlying somatic embryo 
development in spruce have also been investigated (Belmonte 
et al. , 2007; Klimaszewska et al. , 2010). In the current group's 
laboratory, functional characterization has been initiated for 
TFs belonging to a few different classes and has linked two 
MYBs to SCW lignification and related gene deregulation 
(Bomal et al, 2008; Bedon et al, 2010). Spruce overexpres- 
sion lines for both PtMYBl and PtMYB8 displayed ectopic 
lignification and delayed growth and upregulation of genes 
encoding enzymes from the lignin biosynthesis pathway 
(Bomal et al, 2008). 

Depending on research objectives, transient transforma- 
tion can be a valuable alternative offering several advan- 
tages compared with stable transformation. Although it 
cannot reveal morphological phenotypes such as cell-wall 
composition or vascular tissue architecture, it represents a 
rapid approach to obtaining indications of TF-promoter 
interactions. To date, Arabidopsis and tobacco protoplast 
transfections have been extensively used as transient expres- 
sion systems. Nevertheless, protoplast transfection systems 
have been developed for emerging model systems such as 
Physcomitrella patents (Thevenin et al, 2012) that are fast, 
versatile, and robust allowing elucidation of specific TF roles 
in signalling activities and cellular processes. The current 
work aimed to develop a homologous transient transforma- 
tion method in spruce for qualitative testing of promoter-TF 
^raws-activations on a large scale. To achieve this goal, first a 
resource for testing potential interactions of TFs with several 
different promoters in conifers was developed by assembling 
a portfolio of spruce promoters and a large collection of the 
complete coding sequences of members of several families 
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of spruce TFs. This method took advantage of the relative 
uniformity of the cells within embryogenic cell cultures and 
their compatibility with genetic transformation by in vitro 
cultivation with Agrobacterium. The /rara-activation system 
developed here offers a distinct advantage to functional test- 
ing using a heterologous species (such as Arabidopsis) in that 
key cofactors required for TF function are more likely to be 
present in a homologous species. This, however, does not rule 
out the possibility that key factors unique to a specific tissue 
type such as developing secondary xylem may not be present 
in somatic embryonal cells. Nevertheless, it stands to reason 
that a homologous system would better reflect reality com- 
pared to the use of a yeast one-hybrid assay for instance. 

This approach revealed positive interactions of promoters 
resulting from TF fraws-activation. Most of the promoters 
that were tested did not show background expression that 
could result from endogenous expression of TFs. Some TF 
families have been shown to form complexes with other pro- 
teins in order to activate properly their target genes. This was 
found to be particularly important in the regulation of spe- 
cific genes involved in the anthocyanin biosynthetic pathway 
where a number of MYB TFs were shown to form ternary 
complexes with basic helix-loop-helix and WDR type pro- 
teins (Hichri et al, 2011). It is possible that the lack of pro- 
moter-TF ?ra«s-activation is due to the absence of required 
cofactors. However, the method described here could easily 
accommodate for the introduction of a third gene construct 
encoding a specific cofactor, in addition to the TF-coding 
region and the promoter-reporter constructs. Also, any arti- 
ficial /ra«^-activation system is prone to false positives and 
this includes the one described here; the level of transcript 
accumulation of the tested TF is likely to be higher in the 
transiently transformed cells since strong promoters are com- 
monly used. Also, this work cannot exclude the possibility 
that indirect or nonspecific promoter induction could result 
from high accumulation of the TF in the cell. However, this 
study group is quite confident about the robustness of this 
method as the data obtained confirms existing knowledge 
in the literature. As far as is known, the scope of the study 
greatly extends the number of TFs and promoter interactions 
that have been studied in any vascular plant species outside 
of flowering plants. 

Conifer class-IIB NAC-domain proteins display general 
features similar to Arabidopsis 

In Arabidopsis, NAC-domain proteins (SND1, NSTs, and 
VNDs) orchestrate a transcriptional cascade that involves 
R2R3-MYBs and regulates vascular differentiation and SCW 
formation. Functional orthologues of members of this net- 
work described have been identified in poplar and eucalyptus 
(Zhong and Ye, 2009; Zhong et al. , 2010a). For example, pop- 
lar and eucalyptus NAC-domain proteins similar to AtNSTs/ 
AtVNDs have been characterized through their expression 
patterns and restoration of WT phenotype in Arabidopsis 
mutants. Moreover, WNDs and MYB proteins in poplar and 
eucalyptus were shown to be directly involved in the tran- 
scriptional cascade described in Arabidopsis, suggesting a 



relative conservation of a core regulatory network related to 
SCW biosynthesis (Zhong and Ye, 2009; Zhong et al, 2011). 
It was suggested that a similar cascade may also regulate 
SCW formation in conifers (Bomal et al. , 2008), but until the 
current report, functional evidence had only been obtained 
for MYBs (MYB1 and MYB8 from pine and spruce). 

This report shows that PgNAC-4 and PgNAC-7 are the 
only two NAC-domain genes that clustered within the IIB 
subfamily. As such, they represent candidate NAC genes that 
could participate in a cascade regulating vascular develop- 
ment analogous to that described in angiosperm plants and 
trees. The data show that the PgNAC-7 amino acid sequence 
uniquely clusters with AtVND 4/5/6 and with poplar ortho- 
logues PtrVDN 3(A,B) and 4(A,B) and is distinct from other 
VNDs from Arabidopsis and poplar (Fig. 4). Its expression 
is strongly preferential to secondary tissue in both stems and 
roots (Fig. 5) and its transient expression activates the expres- 
sion from promoters of MYBs (MYB1 and MYB8) that have 
been linked to SCW formation, as well as downstream pro- 
teins and enzymes involved in cell-wall assembly (PgCAD, 
Pg4CL, PgTUA-1, PgXTH8-l; Fig. 3). 

In contrast, PgNAC-4 specifically clustered with the 
SMB/BRN subclass genes AtSMB, AtBRNl, and AtBRN2 
involved in cell-type specification and SCW modification 
during root cap maturation in Arabidopsis (Bennett et al, 
2010). PgNAC-4 transcripts were found in all of the Picea 
glauca tissues tested but were clearly preferential (8-10-fold 
difference) to root tips (Fig. 5). In the ?ra«i-activation assays, 
PgNAC-4 gave positive interactions with promoters from 
genes related to SCW formation {PgCAD, Pg4CL, PgDHS2) 
as did PgNAC-7. However, the TFs activated by PgNAC-4 
are different from PgNAC-7. PgNAC-4 was able to positively 
regulate PgHB4 and PgLIM-1 promoters but no significant 
inductions were observed for the PgMYBl and PgMYB8 
promoters (as observed for PgNAC-7). Tissue profiling also 
showed that PgNAC-4, PgCAD, Pg4CL, PgDHS2, PgHB4, 
and PgLIM-1 do not share similar tissue specificity (Fig. 1). 
As only PgNAC-4 transcripts were found in root tips (not 
PgNAC-7), it may be that regulation of PgCAD, Pg4CL, 
PgDHS2, PgHB4, and PgLIM-1 by PgNAC-4 is restricted to 
a small part of the plant. Furthermore, conifer NAC-domain 
TFs from the IIB family may share the same generic ability to 
activate genes related to the SCW as observed in Arabidopsis 
for AtSMB, AtBRNl, and AtBRN2 (Bennett et al, 2010). 

White spruce vascular NAC-domain TF master 
switches 

This work identified two of the eight putative NAC-domain 
proteins as members of the IIB subfamily, whose members 
are known to regulate SCW biosynthesis. Rigault et al. (201 1) 
showed that the family of NAC-domain proteins is one of 
the most underrepresented in Picea glauca, with 36 putative 
members compared with 113 and 162 in Arabidopsis and 
poplar, respectively. This observation is also supported by 
Lang et al. (2010) who showed that TF-encoding genes have 
diversified concomitantly with morphological complexity for 
the development of flowering plants. In angiosperms, water 
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transport is carried out primarily by tracheary elements, and 
mechanical support is provided by vascular and intervascular 
fibre cells. In contrast, both of these functions are carried out 
by tracheids in gymnosperms (Raven et al , 2005). Studies of 
SCW regulation in angiosperms showed that they have tran- 
scriptional regulators that may be specific to particular cells, 
such as NSTs in fibres and VNDs in vessels (Zhong et al, 
2010a). It was suggested that this specificity may be the result 
of the diversification and specialization of vascular tissues in 
angiosperms. Strikingly, the only NAC-domain protein iden- 
tified to date in Picea glauca that could play a master regula- 
tory role (PgNAC-7) is clustered in the VND-like subclade, 
which is specific to vessels. Furthermore, no NST-like gene 
sequences have been identified to date. 

The data presented in this report are helping to delineate 
putative transcriptional networks for the regulation of vas- 
cular tissue differentiation in Picea glauca based on sequence, 
expression, and functional similarity with angiosperm genes 
(Fig. 6). It was previously observed that PtMYB8 and, to a 
lesser extent, PtMYBl may positively regulate genes involved 
in lignification such as Pg4CL and PgDHS2 (Bomal et al, 
2008). PtMYB8 and PgMYB8 are close orthologues of 
AtMYB46 (Bedon et al, 2007), which is a key regulator 
in xylem differentiation that can be directly regulated by 



AtSNDl/NST3, NST1, NST2, VND6, and VND7 (Zhong 
et al , 2007). More specifically, AtMYB46 has been identified 
as a MYB master switch for secondary growth. MYB1 from 
both spruce and pine are highly similar to AtMYB85 and 
AtMYB20 (Bedon et al, 2007) and can be trans-activated 
by PgNAC-7. Moreover MYB20 and MYB85 in Arabidopsis 
have been shown to be positively regulated by AtSNDl/ 
AtNSTl, and AtMYB85 is able to induce the At4CLl pro- 
moter (Zhong et al, 2008). 

An orthologous cascade of ^raws-activations has been 
recently identified in poplar with the characterization of 
WNDs acting as transcriptional master switches (McCarthy 
et al, 2010; Zhong et al, 2010b, 2011; Ohtani et al, 2011). 
Despite the fact that some downstream activated TFs were 
identified as being involved in wood formation, about 13 new 
TFs (not related to secondary growth) have been identified in 
transgenic poplar overexpressing WND genes. They include 
PtrWRKYH and PtrWRKYB that are specifically expressed 
in vascular tissue (Zhong et al, 2011). The current work 
showed that PgWRKY-1 is coregulated with PgNAC-7 and is 
highly preferential to secondary xylem (Fig. 1). Interestingly, 
the frans-activation assays showed that PgWRKY-1, 
whose closest poplar homologues are PtrWRKY12 and 
PtrWRKY13, positively regulate PgMYB8 and most likely 
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PgMYBl (Fig. 3). However, PgWRKY-1 /raws-activation by 
PgNAC-7 and the putative role of this conifer WRKY gene 
in secondary growth remain to be demonstrated. 

Although relationships of orthology between conifers and 
Arabidopsis must be viewed as tentative given that the two lin- 
eages have been separated for around 300 million years, these 
results suggest that PgNAC-7 acts as a first level of transcrip- 
tional control in SCW biosynthesis. Seeing that PtMYB8 
and PgMYB8 can positively regulate PgMYBl, they are 
proposed to function at a second level, upstream of other 
MYBs in a manner analogous to AtMYB46 (Zhong and 
Ye, 2009; Zhao and Dixon, 2011), as depicted in this work's 
model of the transcriptional network regulating SCW bio- 
synthesis and lignification (Fig. 6). PgMYBl is functionally 
related to AtMYB85 with a predominant role in SCW for- 
mation. Recently, a negative -feedback loop has been shown 
for the regulation of AtSNDl promoter activity. Three MYB 
proteins from the Sg4 clade were shown to be downstream 
targets of AtSDNl and were also able to directly repress its 
promoter through binding to a specific cw-motif, thus creat- 
ing a negative-feedback loop (Wang et al. , 201 1). The current 
coexpression study showed that PgMYB3 and PgMYB16 are 
coregulated with genes identified in the putative regulatory 
cascade leading to SCW biosynthesis (Fig. 1). Interestingly, 
these two genes also belong to the Sg4 clade of MYB proteins 
(Bedon et al., 2010). Moreover, the current work showed that 
PtMYB3, the pine orthologue for PgMYB3, is able to repress 
PgDHS2 promoter in the ?ra«s-activation assay, indicating its 
potential repression activity (Fig. 3). 

In summary, PgNAC-7 is cast as a strong candidate among 
the Picea glauca NAC-domain TFs for the role of master 
switch in regulating secondary vascular growth. Overall, 
this work has developed an approach to rapidly delineate 
members of TF families that are potentially linked to physi- 
ological processes through their activity on promoter target 
sequences. This approach is anticipated to accelerate decision 
making to carry out functional analyses such as stable trans- 
formation in Arabidopsis or spruce and more detailed analy- 
ses of the DNA-binding region for a specific TF, as recently 
reported by Bomal et al. (2014). 
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